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Abstract 

Visualization of turbulent flows is a powerful tool to help understand the turbulence dynamics 
and induced transport. However, it does not provide a quantitative description of the observed 
structures. In this paper, an approach to characterize quantitatively the topology of the flows is 
given. The technique, which can be applied to any type of turbulence dynamics, is illustrated 
through the example of resistive ballooning instabilities. 
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I. INTRODUCTION 



Visualization of flow structures plays an important role in providing understanding of the 
turbulence dynamics and some of the properties of the transport mechanisms induced by the 
same turbulence. In the case of turbulent plasmas in toroidal geometry, flow structures may 
be quite complicated [Ij. Eddies form toroidal knots, which through ballooning effects can 
merge in the outer regions of the torus, forming what has been called streamers [2j. In the 
inner toroidal region, there is a fllamentation of these eddies and the resulting flow fllaments 
twist around the torus following the magnetic fleld lines. When these ballooning structures 
remain in the plasma in quasi-steady state, they can induce pseudo-chaotic behavior and 
anomalous diffusion of tracer particles both in the radial and poloidal directions [3l H] . 

Although visualization is a powerful tool in interpreting numerical results, it does not 
provide by itself a quantitative characterization of the flow structures. Here, we take a step 
towards this quantitative characterization by introducing diagnostics that measure some 
topological properties of the structures. The introduction of the proposed techniques will 
be performed through their application to a concrete physical problem: the study of the 
topology of the flows in resistive pressure-gradient-driven turbulence in toroidal plasmas. 
We use the set of reduced magnetohydrodynamic (MHD) equations presented in Refs. [5^ 6j, 
so that the flow is determined by the velocity stream function $ (p, 0, () where (p, 0, () are 
toroidal coordinates. The velocity fleld, V, is given in terms of the stream function as 
V = V$ X b, where b = B/ |B| is a unit vector in the direction of the magnetic fleld. A 2-D 
plot of the constant $ surfaces from a low-/? resistive pressure-gradient-driven turbulence 
calculation is given in Fig. [!} showing the aforementioned streamers in the outer region 
and the fllamentation in the inner region. Algebraic topology can be used to measure the 
complexity of these structures. 

For structures in there are three relevant topological invariants (i.e. quantities which 
do not change under continuous deformations of the structure) : the flrst three Betti numbers 
[7] that we denote by 6o, bi and 62- &o is the number of arc-connected components of the 
structure, bi is the number of (independent) non-contractible loops, and 62 is the number 
of voids. For instance, the solid sphere has 60 = 1, 61 = and 62 = 0, because it has one 
connected component, no voids and any loop can be contracted to a point. However, the 
sphere has 60 = 1, 61 = and 62 = 1, because it has a void. The toroidal surface has two 



2 



independent non-trivial loops (see Fig. [2|) and one void, so that 60 = 1, &i = 2 and 62 = 1, 
but for a solid torus the Betti numbers are bo = 1^ bi = 1 and 62 = because it only has 
one non-contractible loop. 

A crucial point is that the Betti numbers are the rank of commutative groups obtained 
from a succession of morphisms between free commutative groups; they make up the chain 
complex of the space. Furthermore, the Betti numbers also are the dimension of vector 
spaces obtained from a succession of linear maps. As a result, linear algebra methods are 
used to compute them. 

For the complicated flow structures obtained in numerical calculations the determination 
of the Betti numbers is far from trivial. The classical algorithms for the computation of 
the homology of a set [7j (the Betti numbers are the ranks of the homology groups) yield 
exponential run time. Many efforts have been done to reduce the run time and some rather 
efl&cient software has been developed. Here, we use the software package provided by the 
Computational Homology Project (CHomP) [HI [9] in order to diagnose the numerical results 
from 3-D plasma turbulence calculations. The CHomP project has developed many numer- 
ical tools. A description of such tools and some examples of applications can be found in 
the CHomP website ^. 

There are several ways to extract the information needed to construct a chain complex 
of a set in M^. CHomP uses a decomposition of the set in cubes, i.e. it uses a union of cubes 
that approximates our set as flnely as possible. The input data of CHomP are either text 
files listing all those cubes, or a bitmap with each pixel representing a cube. Ref. [9] gives a 
systematic computational approach to the homology of cubical sets, which is the approach 
that we follow here. 

One of the most serious issues that we had to face when calculating the Betti numbers 
with the CHomP software was to determine the right resolution of the 3-D structures, since 
we easily reach the capacity limit of the software. To understand the problems involved in 
finding the proper resolution, we have first investigated, with appropriate analytical models, 
the most eflicient way of providing data about the 3-D structures to the CHomP software 



in order to have converged results for the Betti numbers. This is discussed in Section H 



In Section HI, we present an alternative, somehow indirect method for deducing the Betti 



^ http://chomp.rutgers.edu/ 
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numbers based on 2-D computations with CHomP. In this way we can test the numerical 
convergence and have a verification of the results of Section [11} For simple situations of 



ballooning structures, we give in Section [IV] an analytical calculation of the Betti numbers. 
This provides a way to validate the numerical approach. 

In this paper, we limit the analysis and calculations to a single toroidal ballooning struc- 
ture obtained either through an analytical parameterization or a numerical calculation. 
Studies are under way to calculate the Betti numbers during the time evolution of the 
resistive pressure-gradient-drive turbulence for different beta values. 



II. OPTIMIZATION OF THE CUBE-COVERING OF THE FLOW 

As mentioned in the Introduction, we present our approach to the quantitative study of 
flow topology by means of the particular example of pressure-gradient-driven turbulence in 
toroidal plasmas. The numerical results are the steady state solutions given by the 3D code 
FAR [To], which solves the set of reduced MHD equations j5l|6] as an initial value problem. 
In this code, all fields are Fourier expanded in the poloidal and toroidal angles. We denote 
by m (N) the poloidal (toroidal) mode number. A finite difference representation is used for 
the radial coordinate p. The numerical scheme is implicit for the linear terms, and explicit 
for the nonlinear ones. The evolution is implemented through two half-time steps to ensure 
second-order accuracy in time. Details of the numerical calculations are given in Ref. [llj. 

For very low-/? plasmas, a single toroidal mode, A^, dominates the spectrum jTT]. Here, we 
analyze the topological structure of some of these solutions. Since the functions are defined 
on a grid in p, we either interpolate the values using splines or use a parameterization based 
on the linear structure of these modes. We also use this parameterization as a simplified 
model for these structures and in the verification and validation calculations. 

For a single toroidal mode A", we use a parametrization based on the linear eigenmode 
structure of the resistive ballooning [12j: 
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Here, rriM is the dominant poloidal mode and is the radial location of the singular 
surface with q{pm) = tti/N ^ where q is the safety factor. For a parabolic g-profile, q{p) = 
g(0) + [g(l)-g(0)]p^ we have 



_ m/N-q{0) 

"•"-i <,{!)- <,m- 

Each of the components of the baUooning mode, Eq. ([T]), is characterized by a width Wm- 
When we use the parameterization ([T]) as a simple model in the covering optimization tests 
(see below), we usually take A = 1^ Wb oc, and W . In this way, for a given value 

of A^, the ballooning structure is determined by the single parameter W . When we take the 
FAR numerical results as input for the CHomP software, we either fit the parameters A, 
Wb-j and Wm to the numerical data or use splines to interpolate between radial points. 

The CHomP software computes the Betti numbers of spaces that are union of a finite 
number of n-dimensional cubes, with n < 26; more precisely, spaces that are union of a 
finite number of cubes obtained from of the unit n-cube [0, 1]^ by translation with a vector 
of entire coordinates. We take [0, 1]° to be the origin. Thus the data files used by the CHomP 
program are in the form of a list of cubes, which cover the structure. For a structure in M^, 
the maximum number of cubes that can be used is 2.6 x 10^. Therefore, when we have a 
complicated numerical structure it is wise to find the optimal set of cubes that cover it for 
a fixed number of cubes. 

To optimize the cube selection, the first step is to work on the most efficient coordinate 
system. In the case of magnetically confined plasmas in toroidal geometry, we can use 
Cartesian coordinates in real space, but they are ineflScient because Cartesian cubes do not 
adapt well to the curvature of the toroidal shape. Toroidal coordinates (p, ^, () are better 
suited. Here, p is a radius-like equilibrium fiux surface label, 9 is the poloidal angle and ( 
the toroidal angle. The usage of this coordinate system requires choosing the periodicity 
option for 9 and ( in the CHomP code. When we use toroidal coordinates and for a fixed 
cube size, we need two orders of magnitude less cubes than for Cartesian coordinates with 
the same cube size. These tests were done for a single = 5 structure. 

Most of the numerical tests have been done for toroidal mode numbers in the range 
5 < A^ < 23. The case A^ = 23 is motivated by the results of the nonlinear calculations [llj. 
This mode dominates the toroidal mode spectrum for some of the numerical results. We 
also use lower values of A^ to study the scaling of the numerical properties and because for 
low A^ the calculations are faster and the representation of the structures is simpler. 

The Betti numbers, as any topological invariant, remain unchanged under deformations. 
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Therefore, another way of optimizing the cube selection is stretching or shrinking the coor- 
dinate directions, i.e. making a coordinate transformation (ap, 60,cC), for three convenient 
constants a, b and c. By means of such a change we can choose the number of cubes in each 
direction that cover our torus. In practice, we see that we do not need as much resolution, 
i.e. as many cubes, in the toroidal as in the poloidal and radial directions. It turns out that 
a good choice is Np = No = 2N(^^ where A^^, No and A^^ are the number of cubes in the 
radial, poloidal and toroidal direction respectively. 

Once we have our torus covered with cubes, the next important question is how to 
prescribe that a cube belongs to the structure. For a cube to be in the structure, we require 
that a certain number k of its vertices belongs to the structure; in this case we say that the 
cube is black. 

At a fixed time t, we define a fiow structure as the set of points such that $ (p, 0, C, t) > $0, 
where $o is a constant. A reasonable prescription is that a cube is black if it has k or more 
vertices verifying the condition ^ (p^OX^t) ^ If k is too small there is a problem of 
creation of false loops (overestimation of 61): since the filaments of the ballooning structure 
run quasi-parallel for some toroidal distances and close to each other, if many vertices of 
the cubes lie outside the filaments they can touch and cause the false loops. Fig. [3] shows 
an example of this spurious eflFect in 2-D. In this figure the filaments are represented by 
the regions in red and blue, and we have chosen to consider a cube black if the right-upper 
vertex is either red or blue. We can see the formation of false cycles by the black squares 
(2-D cubes) covering the filaments. 

On the other hand, if k is too large, there are two possible reasons for obtaining false 
results. Firstly, we can have regions in the structure that are thinner than the cubes, causing 
false breaks of the structure, and leading to an overestimation of the number of connected 
components (60) and an underestimation of the number of loops (61). Secondly, the existence 
of filaments very close to each other, with separation of the order of the cube size, produces 
the formation of false loops for reasons analogous to those shown in Fig. [3| In our studies 
we found that this second effect on bi is more important than the first, so we have an overall 
increase of false loops. 

Therefore, in the optimization studies we see that the best selection of k is characterized 
by giving a minimum in bi and the start of an increase in bo. This optimal value of k gives 
the maximum resolution of the structure for a fixed cube size. 
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For the parameters in Ref. [H] and /?o = 0.003, the dominant mode is = 23. We fitted 
the numerical results for this mode using Eq. ([T]) and with this parameterization constructed 
the set of cubes representing the flow structure for k going from 1 to 8. For the flow structure 
with $0 = 0.01 and in Fig.[4| we plotted bo and bi for fixed size cubes, Np = Nq = 2N^ = 200 
as a function of k. The basic picture is as expected. However, the effect of reduction of bi 
is very dramatic as k increases; there is a change of three orders of magnitude in going from 
k = 2 to k = 6. The optimization works for all the different cube sizes that we have tried. In 
Fig. [5] and for the same flow structure as in Fig. [4| we have plotted the results of a systematic 
calculation of bi for different numbers of cubes and for each set of cubes varying the number 
of vertices k used in defining a black cube. In all cases, we have used Np = No = 2A^^. We 
have also plotted in the same graph the value bi = 668 which is the same value for this 



configuration obtained using an alternative approach discussed in Section [IV| We see that 
all values of bi tend to converge to the expected value; however, convergence is very much 
faster for the case of A; = 6 and k = 7 vertices. We have obtained similar results for different 
values of $o- 



III. AN ALTERNATIVE APPROACH TO THE CALCULATION OF THE BETTI 
NUMBERS OF A FLOW STRUCTURE 

Resolution of the numerical calculations is a systematic issue in the determination of 
the Betti numbers. Therefore, it is important to have different ways to calculate them in 
order to do validation of the results. An alternative method to the one discussed in the 
previous section is to consider many toroidal cuts, calculate the 2-D Betti numbers of the 
toroidal sections of the structures and infer the Betti numbers of the 3-D structure. The 
advantage of this approach is that in 2-D calculations we can achieve very high resolution. 
The disadvantage is that it requires performing the 2-D calculations in many toroidal cuts 
if we do not want to miss some of the features of the structure (and it is impossible to know 
a priori how many cuts we need) . 

In 2-D there are only two relevant Betti numbers: the number of connected components, 
6o, and the number of loops, bi. For the structures we are looking for, and in the case of a 
single mode, bi in 2-D is generally zero and all the 2D components are contractible. 

Another important result on Betti numbers is that their alternate sum coincides with the 
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alternate sum of the number of cubes in different dimensions [7] . This number is known as 
the Euler-Poincare characteristic of the set. Assume, for simphcity, that we have a cubical 
space in or M^. By cubes of dimension we mean the vertices of the cubes that form the 
cubical set. By cubes in dimension 1 we mean the edges of these cubes. In the maximal 
cubes have dimension 2 and are squares. In the maximal cubes have dimension 3 and 
their faces are cubes in dimension 2. For example, let A^^, and Ns be the number of 
vertices, edges and squares of the union of a family of squares in M^. Then: 

X = 6o - &i + ^2 = A, - Ae + A,. (3) 

In a ballooning structure, the filaments can be contracted to a 1-D subspace, which is 
just a graph, i.e. a structure formed by points, called vertices, and lines joining them, called 
edges. As this graph is a deformation retract of the filaments, it has the same Betti numbers 
and the same Euler-Poincare characteristic. 

Let us consider for instance a toroidal cut of a given A = 7 structure. In Fig. |6| we 
have plotted in black the regions with $ (p, ^, C = Co) ^ $o for = 0.01 and three values 
of Co- In Fig- 16^5 we can see 10 black connected components. As Co changes the number 
of connected component changes; they break and merge as we move in the C direction. In 
general, given a toroidal cut we will have E well-defined structures. These are identified 
as the first E vertices of our graph. In Fig. [7] we have drawn one such graph with E = 5. 
As we move toroidally, we draw the corresponding number of edges; if we started with E 
vertices we add E edges. Assume that at some toroidal position a filament breaks up into j 
components. This increases the number of vertices and the number of edges by j. Therefore, 
X does not vary. If, on the contrary, at some toroidal position j > 2 filaments merge (Fig.jT]), 
the number of edges increases by j and the number of vertices by 1. So the characteristic 
decreases by j — 1. 

When we complete a full toroidal excursion, the characteristic would decrease each time 
that some filaments merge. Let us assume that there is a decrease in the number of compo- 
nents in going from section s to the next section 5 + 1, i.e. b^^ {s) — bl^ {s + 1) > 0. Then the 
characteristic x will decrease by b^^ (s) — b^^ {s + 1) in going from 5 to 5 + 1. Here, b^^ (s) is 
the 2-D Betti number bo of the structures in the s toroidal cut. Since the system is periodic, 
we find in the last section the E vertices already counted at the beginning. Therefore, the 
characteristic of Euler-Poincare for the full 3-D structure is 
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s 

X = -J2Mi{s), (4) 

where Mi{s) denotes the number of mergers between sections s and 5 + 1, i.e. Mi (s) = 
{^) ~ ('^ + 1) when this number is positive and otherwise. S is the total number of 
toroidal sections and Mi{S) is the number of mergers between sections S and 1. 

In these calculations we are assuming that mergers and bifurcations do not occur si- 
multaneously. This is true for the analytical models that we use provided that we take a 
sufficiently large number of toroidal sections. 

In case of a single ballooning mode and for low values of $0? &o = 1 and 62 = 0, therefore 
bi = 1 + ^5^1 Ml (s). However, this is not true in general and when we use this method in 
more complicated structures, we will only make the determination of the characteristic of 
Euler-Poincare. 

As pointed out above, the main problem with this approach is that, if we do not use 
enough toroidal sections, we might miss some of the filaments mergers. Therefore, this type 
of calculation can give us a lower bound to the Betti numbers. As an example, we calculate 
the bi number for the same = 23 structure as in Fig. [5] using the direct computation 
explained in Section [TT] and the method just described. In Fig. |8| we compare the results 
for both methods. We have plotted bi calculated with the method of Section [TT] in the case 
of highest resolution, Np = No = 2N(^ = 1200, as a function of $0 for a range of values in 
which there is not much variation on the results. We have also plotted the results given 
by the method discussed in this section for a 2-D Cartesian grid of 8000 x 8000 and using 
different numbers of toroidal cuts. We can see a relatively good agreement between the 
two calculations, although a very large number of toroidal cuts are needed to get converged 
results. 

For a single toroidal mode the approach of this section can work very well. However, when 
we have a broad spectrum of toroidal modes the situation is more complicated, essentially 
due to the fact that we can have loops in the 2-D cuts. The break of a 2-D loop does 
not change the Euler-Poincare characteristic. The merger of two portions of a connected 
structure to form a new loop decreases x by 1. In case that a loop contract to a point, the 
characteristic increases by 1, but this phenomenon has never been observed in the structures 
we are studying and we are not going to take it into account. Therefore now we should 
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calculate the characteristic as follows. Let us assume that we consider S toroidal cuts and 
denote by (s) and b^^ (s) the Betti numbers of the 2-D structures in the s toroidal cut. 
We define 



Mi{s)= { V 7 I . u V . u V . 



if 
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X = -J][Mi(s) + M2(s)]. (7) 

This is the expression we use for the more general nonhnear flow structures that we flnd 
in the 3-D nonhnear calculations. 



IV. AN ANALYTICAL MODEL FOR THE BETTI NUMBERS OF THE BAL- 
LOONING MODES 

Let us calculate the Betti numbers for a ballooning mode with the structure given by 
Eq. ([T]). We will use the approach described in the previous section to do the calculation of 
the Betti numbers. 

One way of understanding the structures in 2D given by Eq. ([T]) for constant values of ( 
is by flrst making a model for the different components of the ballooning mode. We identify 
each component m with the position of its m maxima. They are located at the rational 
surface m/N and at the poloidal positions 

^ Ak + ln m,m — 1 

m 2 m 2 2 

Let us consider a fleld period, tt/ (2A^) < ( < Stt/ {2N). We begin the fleld period at 
tt/ (2A^) because at this toroidal cut the distribution of structures has maximum symmetry. 

For simplicity we assume that $ has been normalized in a way that its maximum value 
is 1. We consider a range of values for $o ranging from 0.001 to about 0.01. In the next 
section, we will discuss in more detail the range of validity. The regions of local maxima of $ 
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tend to merge forming streamers. This merger can be realized by joining the local maxima 
at different radii without crossing lines of minima. For N = 7^ and in Fig. [9} we show the 
resulting structures from the merger of the maxima. Note that the lines joining the maxima 
correspond to the black structures in Fig. [6| 

In Eq. g, mo and rui are the lower and upper hmits of the sum over poloidal modes. 
Assuming that mo is even and mi is odd and using the model for the structures represented 
in Fig. |9| we count the number of structures in this toroidal cut. For very low $0, we have 
the following structures: 

1) The one at ^ = 0, which is common to all. 

2) The (mo — 2)/2 pairs before the merge of the positive and negative branches. 

3) The (mi — mo + l)/2 structures that merge across the x-axis. 



The total is 



Ns. = !!!i±|!illi. (9) 



In other toroidal cuts, the number of structures varies due to bifurcations and mergers 
between structures. Let us look closely at these phenomena for that simple model. 

The process of bifurcation and merger of the structures seems to follow a general pattern. 
All the (mi — mo + l)/2 structures that result from mergers across the x-axis on the left-hand- 
side bifurcate after a toroidal rotation of 7t/{2N). Therefore, the number of components at 
this toroidal angle reaches its maximum value mi. After another rotation of 7r/(2A^), we 
see that some of the bifurcated structures merge. In merging, they do not join the same 
structure as before, but the next one. This allows the rotation of the whole system. In 
the second half of the field period, the process reverses. After a field period is completed 
the structures are the same as the initial ones, but all maxima have performed a poloidal 
rotation of 27r/m. 

In short, the number of structures in each quarter of a field period is 

mo + mi-1 mo + mi + 1 mo + mi - 1 , 
^ >mi^ >mi^ . (10) 

All bifurcations happen on the negative x-axis and so do mergers. The reason for this 
behavior is that the filaments in the structure move from left to right in the upper half plane 
and from the right in the lower half plane. Once a structure crosses the axis it is torn apart. 
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Here we have done the discussion of the structures for given parity of mo and mi. For the 
general case, results are given in Table |T| 

For a low range of values of $o and for the parameterization of Eq. ([T]) , there is only one 
connected component in 3D and no voids. Therefore, we can use Eq. ^ to get bi. Namely, 

60 = 1, 61 = 1 + A^(mi -mo), and 62 = 0. (11) 

The different parity combinations listed in Table |T] give the same value for the bi Betti 
number. 

This (approximate) analytical calculation of the Betti numbers is very useful for the 
validation of numerical calculations. An example is shown in Fig.flOlfor an = 5 structure. 



We can see the numerical results converging to the value given by Eq. (11) as the number 
of cubes in the radial direction increases. 

For the numerical calculations in which the A = 23 mode dominates, the relevant of 
poloidal mode numbers (modes such that > $0) is from mo = 24 to mi = 53. Using 



Eq. (11), we obtain bi = 668, this value is consistent with the numerical results from the 



CHomP code discussed in Sect. HIl 



V. SUMMARY AND CONCLUSIONS 



The results from Sections [TT] and [TTT] show that we have a double approach in calculating 



the Betti numbers of a turbulent flow structure. One is by constructing a 3-D optimized 
covering of the structure with cubes. As the size of the cubes decreases and the number of 
cubes increases, the Betti number computations converge. 

Consider the A = 23 structure of Ref. ^Ij. The Betti number bi computed using the 



method of Section II is shown in Fig. 11 for different values of $o/$Max- Here, $Max is 
the maximum value of $ in the whole plasma volume. As the number of cubes increases 
bi converges from above. In these calculations, A^ is a multiple of A and A^ = Nq are 
multiple of 200 and close to 2A^. The same calculation performed along the lines of Section 



III involves the determination of bi from a number of 2-D toroidal cuts through Eqs. ^ and 



([6]). The results are shown in Fig. 12, Here, Np = No = 2400, and the number of toroidal 
cuts is always a multiple of A. As the number of toroidal cuts increases, bi converges 



from below as it is shown in Fig. 12, The only exception is the case with A^ = 2392 and 
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$o/$Max = 0.002, that could indicate that we need higher resolution in the (p, ^) plane for 
that value of A^^. This convergence behavior gives a nice way to estimate this Betti number 
in more complicated structures. By applying the two techniques for calculating fei, we can 
have an upper and a lower bound to its real value. 

The relevant range of values of $o is another important consequence of these studies. It 
seems that the optimal range is < $o/$Max ^ 0.01. The Betti numbers are practically 
constant in this range. In Fig. 13, we show the Betti number bi as a function of $o/$Max 



for three different values of A^. Note that bi is constant and has its maximum value in this 
low range of values of $o- The figure also shows that the value of bi in this range of $o is 



the one given by the analytical calculation of Section |IV[ For all these reasons, the choice 
of this range for $o seems to be the best to characterize the structure. As $o increases, the 
fiow filaments become narrower and break, causing a drop of bi. 

In conclusion, we have shown that it is possible to give a quantitative characterization of 
the turbulent plasma fiows by determining the Betti number corresponding to structures of 
constant values of the stream function. The CHomP project offers a broad set of tools to 
carry out the numerical calculations. 

Having a good resolution of the surfaces is essential. The intricate topological structures 
of the ballooning modes can lead to many false values for the number of loops and careful 
convergence studies are required. 

We have used two approaches in those calculations: 

1) Constructing a full 3-D file of cubes describing the structure. We have found a way of 
getting optimal resolution for the determination of the cubes. This approach tends to 
give an upper bound of the 3-D bi Betti number. 

2) Using multiple toroidal cuts and determining the Euler Poincare characteristic from 
the number of mergers of 2-D components. This approach gives a lower bound to the 
3-D bi Betti number. 

For a simple fixed A^ ballooning modes, the two approaches give in general the same value. 
For more complicated structures, we will use the two approaches in order to get upper and 
lower bounds for the Betti numbers. 

The next step in this research is to characterize the variation, in time and (3^ of the 
Betti numbers for the fiows obtained in 3-D resistive pressure-gradient-driven turbulence 
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calculations. 
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TABLE I: Number of 2-D connected components at different toroidal sections. 



tt/{2N) tt/N 37r/(2iV) 27r/N 57r/{2N) 

mo even, mi even mi±nii ^1 mo±nii ^1 mo±mi 

mo odd, mi odd ^^^^ mi ^^^^ mi ^^^2±^ 

mo even, mi odd ^O+^l-l rno+mi + l mo+n^i-l 

mo odd, mi even iM+m+i mo+mi^ mo+mi+i 
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FIG. 1: Contours of $ = Constant on the toroidal cross section with C = 0. Streamers form in the 
outer region (right) of the torus. 
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FIG. 3: 2-D representation of the covering by squares of two filaments running nearly parallels. 
The regions in red and blue represent the filaments. The covering cubes are black if the right-upper 
vertex is either red or blue. We can see the formation of false cycles by the black squares (2-D 
cubes) covering the filaments. 



18 




number of vertices k 

FIG. 4: Betti numbers 60 and bi for an = 23 flow structure as a function of fc, the number of 
vertices in the structure that defines a black cube. The calculation is for a fixed size cubes with 
Np = Ne = 2Nc^ = 200. 
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FIG. 5: Betti number hi for the same = 23 flow structure as in Figure [4] as a function of A:, the 
number of vertices in the structure that define a black cube. The calculation shows the convergence 
as the number of cubes increase for cubes coverings verifying Np — Nq — 2Nq. 
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FIG. 6: Regions with $ (p, 6>, C = Co) > $o for a TV = 7 flow structure. Here $o 

plotted three values of Co- a) Co = 7r/14, b) Co = and c) Co = 37r/14. 
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FIG. 7: Graph representing the changes in the connected components as we move in the toroidal 
direction. The nodes represent the connected components in a given toroidal cut. 
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FIG. 8: Comparison of the results for the Betti number bi as a function of $o calculated with the 
methods of Section [IT] and Section IIL For the first, using Np = Nq = 2N(^ = 1200, and for the 
second using a 2-D Cartesian grid of 8000 x 8000. The latter has been done for different number 
of toroidal cuts. 
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FIG. 9: Structures from the merger of the maxima of $ at different radii for the same N = 7 flow 
structure as Figure |6] for the same three toroidal cuts shown in that figure. Note that the lines 
joining the maxima correspond to the black structures in Figure |6| 
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FIG. 10: Numerical results of the Betti number 61 for a = 5 flow structure converging to the 



value of the analytical model of Section IV, Eq. (11 ), as the number of cubes in the radial direction 
increases. 
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FIG. 11: The Betti number bi computed using the method of Section [Tl] for a numericahy calculated 
= 23 flow structure [11 J and different values of $o/^Max- As the number of cubes increases bi 
converges from above. The horizontal line corresponds to the converged value bi = 668. 
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FIG. 12: The Betti number bi computed using the method of Section [III| for the same cases as Fi^ 
11 The number of toroidal cuts is always a multiple of = 23, and Np = Nq = 2400. 
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FIG. 13: The Betti number bi as a function of ^>o/^Max for three different values of N. It shows 
the range of values of $o in which bi is constant. The figure also shows that the value of bi in this 



range of $o is the one given by the analytical calculation of Section IV 
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